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auxiliary  variable,  usually  representing  another  physically  important  quantity, 
and  writes  the  differential  equation  as  a  lower  order  system.  One  then  con¬ 
siders  Ritz-Galerkin  approximation  schemes  based  on  a  variational  formulation 
of  this  lower  order  system,  thereby  obtaining  direct  approximations  to  both 
the  original  and  auxiliary  variables.  Three  particular  mixed  methods  for  the 
approximate  solution  of  the  biharmonic  problem  are  examined  in  detail. 
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ANALYSIS  OF  MIXED  METHODS  USING 
MESH  DEPENDENT  NORMS 

I.  BabuSka  J.  Osborn^,  and  J.  PitkSranta  ^ 


1.  Introduction 

In  15]  Brezzi  studied  Ritz-Galerkin  approximation  of  saddle  point  problems  arising 
in  connection  with  Lagrange  multipliers.  These  problems  have  the  form: 

Given  f  €  V'  and  g  c  II’’  ,  find  €  V  »  W  satisfying 

a(u,v)  +  b(v,i^)  »  (f,v)  V  V  f  1/ 

(1.1) 

b(u,>^)  ^  W  ^  i  W  , 

where  V  and  it’  are  real  Hilbert  spaces  and  a(*,*)  and  b(',*)  are  bounded 
bilinear  forms  on  K  x  (/  and  V  x  W  ,  respectively. 

Given  finite  dimensional  spaces  V,  c  (/  and  w.  c  u;,  indexed  by  the  parameter  0  ^  h 

h  h 

the  Ritz-GalerJcin  approximation  (u.  ,<i.  /  to  (u.ip)  is  defined  as  the  solution  of  the 

h  n 

problem: 

Find  (u.  ,4i.  )  e  V.  *  W.  satisfying 
h  h  h  h 

r  a(u,  ,  v)  +  b(v,\(i  )  »  (f,v)  V  V  €  V 

(1.2)  [  h  h  h 

I  b(Uj^,^)  »  V  e  . 

The  major  assumptions  in  Brezzi 's  results  are 

(1.3)  sup  ^  YqIIuIIj,  V  u  £  and  v  h  , 

n 

where  Yq  >  0  is  independent  of  h  ,  and  Zj^  *  {v  £  Vj^:b(v,i^)  =0  V  •f  c 

(1.4)  sup  >  kJv’ll,.  V  £  W  and  V  h  , 

veV^  ~ 

n 

where  )£q  >  0  is  independent  of  h  .  Using  (1.3)  and  (1.4)  Brezzi  proves  the  following 

error  estimate  for  the  approximation  method  determined  by  (1.2): 


(l.S)  Ku-u.ll|.  +  ll’h'I'JlM  T  Cdnfllu-yllj.  +  infllij:-nll^)  V  h, 

xeVh  n£Wj^ 

where  C  is  independent  of  h  . 
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In  (1,21  BabuSka  studied  Ritz-Galerkin  approxiraation  of  general,  variationally 
posed  problems.  The  main  result  of  [1,2],  as  applied  to  (1.1)  and  (1.2),  is  that  (1.5)  holds 


provided 


(1.6) 


where  >  0  is  independent  of  h  .  It  is  clear  from  [1,2]  that  (1.3)  and  (1.4)  hold  if 
and  only  if  (1.6)  holds.  (1.3)-(1.4)  or,  equivalently,  (1.6)  is  referred  to  as  the  sta¬ 
bility  condition  for  this  approximation  method. 

The  results  of  [1,2,5]  can  be  viewed  as  a  stratagy  for  analyzing  such  approximation 
methods:  the  approximation  method  is  characterized  by  certain  bilinear  forms,  norms  (spaces), 
and  families  of  finite  dimensional  approximation  spaces,  and  if  the  method  can  be  shown  to 
be  stable  with  respect  to  the  chosen  norms,  then  the  error  estimates  in  these  norms  follow 
directly,  provided  the  bilinear  forms  are  Ixiunded  and  the  approximation  properties  of 

and  are  )uiown  in  these  norms.  These  results  can  be  used  to  analyze, for  example,  certain 

hybrid  methods  for  the  biharmonic  problem  [5,6].  The  results  of  [1,2]  have  also  l>een  used 


to  analyze  a  variety  of  variationally  posed  problems  that  do  not  have  the  form  (1.1). 

There  are  other  problems  of  a  similar  nature,  however,  where  attempts  at  using  the 
results  of  [1,2,5]  were  not  entirely  successful  since  not  all  of  the  hypotheses  were  satis¬ 
fied:  specifically,  the  Brezzi  condition  (1.3)  or,  equivalently,  the  BabuSIca  condition 
(1.6),  is  not  satisfied  with  the  usual  choice  of  norms,  i.e.,  the  approximation  methods 
for  these  problems  are  not  stable  with  respect  to  the  usual  choice  of  norms.  This  is  the 
case,  for  example,  in  the  analysis  [7l  of  the  Kerrmann-Mlyoshi  [15,16,20]  mixed  method  for 
the  biharmonic  problem.  In  the  analysis  of  this  method  a  natural  choice  for  both  II  *11 
and  ll•ll(^r  is  the  1st  order  Sobolev  norm;  however,  this  method  is  not  stable  with  respect  to 
this  choice  of  norms.  As  a  result  of  this  difficulty  the  error  estimates  obtained  in  [5 1 
are  not  optimal.  A  similar  difficulty  arises  in  the  analysis  of  the  Herrmann-Johnson  [15, 

16,17]  and  Ciarlet-Raviart  [9]  mixed  methods  for  the  biharmonic  problem.  In  later  work  of  Scholz  ' 

1 

[23]  and  Rannacher  [22]  optimal  error  estimates  were  obtained  for  the  mixed  methods  of 
’arlet-Raviart  and  Herrmann-niyoshi.  In  a  forthcoming  paper  Falk-Osljorn  [12]  develop 
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abstract  results  froca  which  optimal  error  estimates  for  these  and  other  problems  can  be 
derived.  However,  in  neither  the  work  of  Scholz  [23],  Rannacher  [22],  or  Falk-osborn 
[12]  is  the  systematic  approach  of  Brezzi  and  BabuSka  used. 

It  it  the  purpose  of  this  paper  to  analyze  mixed  methods  for  the  biharmonic  problem  by  means 
of  the  results  of  Brezzi  and  BabuSka.  This  is  done  by  introducing  a  new  family  of  (mesh 
dependent)  norms  with  respect  to  which  the  above  mentioned  mixed  methods  (Ciarlet-Raviart , 
Herrmann-.v.iyoshl ,  Herrmann-Johnson)  are  stable.  Once  the  stability  condition  has  (seen 
checked  and  the  approximation  properties  of  the  subspaces  Vj^  and  have  (seen  determined 

in  these  new  norms,  the  error  estimates  in  these  norms  follow  immediately  from  the  abstract 
results  of  Brezzi  and  BabuS)ca.  Error  estimates  in  the  more  standard  norms  are  then  obtained 
by  using  the  usual  duality  argument.  The  results  of  this  paper  were  announced  in  [21].  We 
also  note  that  the  methods  employed  in  this  paper  have  been  applied  to  two  point  boundary 
value  problems  in  [3] . 

Section  2  contains  a  reviev  of  the  convergence  results  of  Brezzi  and  Babu5)ca.  In 
Section  3  we  introduce  and  study  the  mesh  dependent  norms  and  spaces  used  in  the  analysis 
in  this  paper.  In  Section  4  we  treat  three  examples  previously  analyzed  in  the  literature 
and  show  how  error  estimates  can  be  derived  from  the  abstract  results  in  Section  2,  used  in 
conjunction  with  the  mesh  dependent  norms  introduced  in  Section  3.  These  exeunples  are  all 
mixed  methods  for  the  biharmonic  problem.  The  error  estimates  in  the  st^mdard  norms  that 
are  obtained  in  the  present  paper  and  those  obtained  in  [12],  using  different  techniques, 
are  the  same. 

Throughout  this  paper  we  will  use  the  Sotxjlev  spaces  h"”  =  H™  (0)  ,  where  (2  is  a 
convex  polygon  in  the  plane  and  m  is  a  nonnegative  integer.  On  these  spaces  we  have  the 


seminorms  and  norms 

|v! 


'm 


^  .1  /  •! 


and 


llvll  =  llvll  «  (  ?  f  iD'Vl^dx)^’^^  . 


a  <  m  n 


denotes  the  subspace  of  ^"((2)  of  functions  vanishing  together  with  their  first 
m-1  normal  derivatives  on  r  «  9(2.  We  also  use  the  spaces  H  "((2)  =  (Hq((2))'  (the  dual 
space  of  H^(ri))  with  the  norm  on  h”'”((2)  taken  to  lae  the  usual  dual  norm. 


2.  Abstract  Convergence  Results 


In  this  section  we  review  certain  results  on  the  approximate  solution  of  saddle 


point  problems 


Let  and  W  be  real  Hilbert  spaces  (indexed  by  the  parameter  h  ,  where  0 


with  norms 


where  K.  and  K,  are  consteuits  that  do  not  depend  on  h 


We  consider  the  following  problem,  referred  to  as  problem  P 


find  (u,i())  £  V 


Given 


)  denotes  the  pairing  between  U  and  its  dual  space  i'' 


We  shall  consider  this  problem  for  a  subclass  of  data,  i.e 


where  D  is 


We  assume  that  P  has  a  unique  solution  for  all  (f,g)  ■:  D 


We  are  interested  in  the  approximate  solution  of  P  .  Toward  this  end  we  suppose  we 


are  given  finite  dimensional  spaces  V  c  (/  and  W  c  (t' 


following  problem,  referred  to  as  problem  P 


Given  (f,g)  e  D  ,  find  (u  ,0  )  e  V.  x  w  satisfying 


We  now  regard  u  as  an  approximation  to  u  and  ij;  as  an  approximation  to 


(2.6) 


V 


and  V  h  , 


where  >  0  is  independent  of  h  .  We  now  state  the  fundamental  estimate  for  the 

errors  u  -  u.  and  fli  ~  ■ 

h  h 

Theorem  1  (Brezzi  (51).  Suppose  (2.1),  (2.2),  (2.5)  and  (2.6)  are  satisfied.  Then 

Problem  has  a  unique  solution  (u.  ,0.  )  for  each  h  and  there  is  a  constant  C  , 

h  h  h 

independent  of  h  ,  such  that 

(2.7)  (lu-uJI,,  +  <  C(  inf  llu-xll,,  +  inf  )  v  h  . 

h  h  ^ 

(2. 5) “(2.6)  is  referred  to  as  the  stea>ility  condition  for  this  approximation  method. 

In  many  applications  of  Theorem  1  the  spaces  V  and  ti  and  the  forms  a,  and  b 

n  h  nr. 

do  not  depend  on  h  ,  i.e.,  I'  «  1/  and  W  =  W  are  fixed  Hilljert  spaces  and  a.  =  a 

h  h  n 

and  b  =  b  are  fixed  bilinear  forms  and  (/  x  (/  and  V  ^  W  .  The  space  V.  and  w. 

h  h  h 

typically  are  spaces  of  piecewise  polynomials  with  respect  to  a  triangulation  of  some 

domain  by  triangles  of  size  less  than  or  equal  to  h  and,  of  course,  depend  on  h  .  In 

the  applications  in  this  paper,  both  the  spaces  (/.  ,  W.  and  v.  ,  W.  depend  on  h,  i.e., 

h  h  h  h 

are  mesh  dependent;  the  constants  Kj^,K2,  y^,  and  k^,  however,  will  be  independent  of  h 
(cf.  [2,  Cp.  7]).  In  these  applications  the  solution  (u,i^)  of  (2.3)  is  independent  of  h 
and  lies  in  x  (</j^  for  all  h  .  Thus  the  estimate  (2.7)  yields  convergence  estimates 
for  families  ^'^h^  satisfy  an  approx imobility 

assumption.  For  typical  finite  element  applications  this  would  involve  the  assumption  that 


inf  llu-xll  |.  and  inf  lli(i-nlL,  tend  to  zero  as  h  tends  to  zero. 
X  e  Vjj  ‘'h  ncWj^  h 


3.  ‘'.esh  Dependent  JJorms  and  Spaces 


In  this  section  we  describe  the  mesh  dependent  norms  and  spaces  we  shall  use  in  the 


be  a  convex  polygon  in  the  plane.  For  0  <  h  <  1  we  let 


angulation  of  '  by  triangles  T  of  diameter  less  than  or  equal  to  h 


family  of  triangulations  {T  )  satisfies  the  minimal  angle  condition,  i 


constant  c  such  that 


max 

T?T. 


where  h  is  the  diameter  of  T  and  c_  is  the  diameter  of  the  largest  circle  contained 


there  is  a  constant  i  >  0  such  that 


and  is  quasi-uniform,  i.e 


define  the  norm 


ST'*'  n  ST  is  an  interior  edge  of  the  triangulation 


where  -r  is  the  unit  normal  to  T'  exterior  to  T 


is  a  boundary  edge  of  T 


and  then  define  H, 


to  be  the  completion  of  H  (Q)  with  respect  to 


be  identified  with  L  (Q)  ®  L, ( 


we  note  that  norms  similar  to  II  *11 


have  been  used  in  a  different 


manner  in  Douglas-Dupont  [11]  and  Thomas  [26) 
For  h  1  a  fixed  integer  we  define 


where  P  is  the  space  of  polynomials  of  degree  k  or  less  in  the  variable  x  and  x 


t  is  clear  that  S  is  contained  in  H 


We  now  prove  several  lemmas  that  are  fundamental  to  the  analysis  of  this  paper 


These  proofs  are  all  closely  related  to  the  ideas  used  in  the  proof  of  the  Bramble-Hilbert 
lemma  N).  Prior  to  stating  the  first  of  these  lemmas  we  describe  the  notation  we  will 


use  and  state  some  well-known  results  that  will  be  used  in  the  proofs 


Let  T  be  an  arbitrary  triangle  and  let  f  be  the  reference  triangle  with  vertices 


(0.0),  (1,0),  and  (0,1).  Then  there  is  an  invertible  affire  mapping  F  (x)  =  B  x  +  b. 


F(x)  =  Bx  +  b  such  that  T  =  F  (f ) .  This  mapping  leads  to  the  correspondence  x  e  T  "  x 


F  (x)  <  T  (between  points  in  T  and  points  in  T  and  the  correspondence  (v:T->-R) 


T  -*■  R)  Ijetween  functions  defined  on  T  and  functions  defined  on  T  .  Note 


that  v(x)  =  v(x) 


It  is  easily  seen  that 


If  v  =  u(x)  denotes  the  outward  unit  normal  to  3T  at  x  and  v  =  v(x)  is  the  outward 


unit  normal  to  3T  and  x 


then 


area  of  T  and  |T!|  denotes  the  length  of  T 


The  seminorms  v 


where  llBlI  is  the  norm  of  B  induced  by  the  Euclidian  vector  norm  (cf.  [8  ,  Theorem 


(of.  [8,  Theorem  3.1.3]).  We  also  note  that  |det  b| 
there  is  a  constant  C  =  C(T)  such  that 


Finally  we  remark  that 


(cf.  [8,  Theorem  3.1.1)) 


There  is  a  constant  C  such  that 


Proof.  It  is  sufficient  to  show  that 


are  both  norms  on  the  finite  dimensiona 


and  hence  there  is  a  constant 


such  that 


Let  T  c  T  and  suppose  T  is  the  image  of  T  under  the  mapping  F (x)  =  Bx  +  b  .  Then 


(3.2), and  (3.6),  we  see  that  for  any  u  €  P  we  .have 


max 


£  C(f)  max|T;||det  b| 


<_  C(T)  max  |t:  I 


Therefore 


h  /  ju  l^ds  £  h  /  (u  I  ds 


‘h 


TeT  3T 
h 


<  C(f)  o  T  y  Hull 
-  „  ‘■t  0« 


T  eZ 


<  C  (T)  0  T  II  ull 


o,n 


for  all  u  £  S, 


Lemma  2 .  There  is  a  constant  C  such  that 


"“"2,h  ^’^"""i.n 


V  U  £  P, 


Proof.  Since  {T.  }  is  quasi-uniform  it  is  well-knovm  that 

""  ”  h 


I 


Tft. 


Iiul|2^^  <Ch-^ 


V  u  ^  P, 


Thus  it  is  sufficient  to  show  that 

h’M  |J  ds  1C  h"2  llul|2 


i,n 


V  U  £  P^  . 


(J^  |ul^dx  +  J[  |v.u|^ds)^^^  and  llullj^  .j  are  both  norms  on  the  finite  dimensional 


3T 

space  and  hence  there  is  a  constamt  C(f)  such  that 


E(u)  5  /  |v.ul  ds  1  C(T)llullj^^;^ 


Clearly 


af 


E(u  +  p)  -  E(u)  V  p  £  Pp 


V  U  f  P,  (T)  . 

k 


Thus 


E(u)  =  E(u  +  p)  1  C(T)  Hu  +  pHj^.j; 
and  hence,  using  (3.9),  we  have 

E(u)  <  C(T)  inf  Hu  +  pH 


PfPr 


1,  r 


P  £  Pr 


<  C(T)  u 


1,T 


Now  let  T  £  T.  and  assume  T  is  the  image  of  T  under  the  mapping  F(x)  =  B.x 
h 

Then,  using  (3.1),  (3.2),  (3.4),  (3.6),  and  (3,8),  we  see  that  for  any  u  e  Pj^  ve  have 


J,  /  I  t(V  u)  (x))''v(x)  l^ds 

^  m  I  ^ 


max 


<  C(T)llB 


<  C(T)|}B 


max 


This 


There  is  a  constant  C  such  that 

inf  Hu-'  II  £  C  h^  |u| 


where  1  <  r  and  1  <  i  <  min{r,k+l) 


Proof.  We  define  two  interpolation  operators  that  will  be  used  in  the  proof.  For 


Then  for  u  £  H  (.1)  we  let  7  u  €  S 


'or  u  £  H'''(r>)  we  define  the  interpolant  in  a  different  manner.  Here  we  consider 


Now,  following  Clement  [10]  we  define  I  u  by 


In  this  case  we  obtain  the  desired 


By  the  standard  approximability  results  for  we 


have 


Thus  it  is  sufficient  to  show  that 


Suppose  u  £  h‘’’(T)  and  set  E(il)  -  /  -  Uul^ds.  By  the  trace  theorem  and 

3T 

the  Sobolev  imbedding  theorem  we  have 

E{u)  <  C(T)  - 

—  ll,T 

and  since  E(u  +  p)  -  E(u)  V  p  £  ,  we  thus  have 

E(u)<_C(T)  inf  llu  +  pllj^  '  i  C(f)  |u|j^  ^  . 

P"Vl 

Now  let  T  £  r  be  the  image  of  T  under  the  mapping  F(x)  -  Bx  +  b  . 
n 

Then 


/|u 

3T 


I  ul^ds  1  ^i  /  I"  - 
T 

^  max|T;  I  /Ju  -  I^ul^ds 


3T 


<max  1t|1  C(T)1u1^^- 


£  max  |t:  lc(T)  |det  bI'^  llBlI^'lujj  ^ 
i  ^ 


4h_ 


<  C(T)  T| 


21 


'i,r 


C'r)  jrl  4  o  |^|2 


^  p,; 


2t 


e,T 
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This  completes  the  proof  for  the  case  r ,  li.  ^  2 


[lOJ  has  shovm  that 


By  a  slight  modification  of  the  proof  in  [91  we  obtain 
ti.  t  I..  T  -  I'  kK.I 


The  desired  result  now  follows 


Lemma  4.  There  is  a  constant  C  such  that 


where  2  <  r  and  2^1  £  min(r,k+l) 


Note  that 


Since,  by  standard  approx imability  results  we  have 
r  11  t  ._ii2  ^  _  .2i-4  I.  I  2 


We  next  observe  that 


Now  let  T  «  T  be  the  image  of  T  under  the  mapping  F(x)  «  Bx  +  b 


Then 


/  (u-I^u)!^ds=  -  I^u)  l^v(x)  !^ds 

3  t; 


y  /  I  (b'-'-)*'  v.(u  -  i-u;  rds 

^  ip »  *  T 

i 


i  li  i(B'S''v.(u  -  I-G)  (i)  |^|t'  Ids 

<  Ib‘^1^  max  It: I  /  |V.(u  -  I-u;  I^ds 
1  ai;-  ^  ‘ 


<  C(f)lB‘^l|2  ™ax|T'||u|2^- 


-  12  rv,  ;  24 


‘'T  rT 


C(T)|T|  4  .3 

- 24“  °  \  l"l4,T 

^  Pt 


Therefore 


3(U-U  )  35 

/  <3s  £  y  / 


TcT.  3T 
h 


4  .24-3  I  1 2 

£  C(T)  o  h  |u|^^^  , 


which  completes  the  proof. 


4.  Applications 


In  this  section  we  analyze  three  mixed  methods 


a)  Ciarlet-Raviart  method 


Consider  the  biharmonic  problem 


is  a  convex  polygon  in  the  plane  and  g  is  a  given  function.  If  g  £  H  (fJ) 


then  there  is  a  unique  solution  e  H  (fi)  of  (4.1).  In  addition  the  following  regularity 


result  is  )(.nown  for  this  problem:  If  g  c  H  (P.)  ,  then  i);  £  H  (P)  n  H  (P)  and  there  is 


a  constant  C  such  that 


Using  the  well-)cnown  correspondence  between  the  biharmonic  problem  and  the  Stokes  problem 
this  regularity  result  can  be  deduced  from  the  regularity  result  for  the  stokes  problem 


proved  in  [  18]  .  We  assume  g  £  H  (P)  throughout  this  section 


We  now  seek  an  approximation  to  the  solution  il>  of  (4.1)  toy  a  mixed  method,  i.e., 


we  introduce  an  auxiliary  variable  (u  s  -Aii  for  the  method  of  this  subsection)  ,  write 


(4.1)  as  a  second  order  system,  cast  this  system  into  variational  form,  and  then  consider 


the  Ritz-Galerkin  method  corresponding  to  this  variational  formulation 


The  desired  variational  formulation  of  (4.3)  is  obtained  by  multiplying  the  1 —  equation 


integrating  the  resulting 


the  2 —  equation  toy  v  c  H. 


and  integrating  the  first  one  by  parts  over  each  T  £  T  .  By  means 


of  this  process  we  arrive  at  the  following  problem 


n  H  )  satisfying 


find 


I  u  V  dx  -  J  I  V  i.  i  dx  -  j  v(J  )ds  =  0  V  v 

n  T'T  t  r, 

n  n 

(4.4) 

y  /  u  A  ^  dx  -  /  u  (j|^  )  ds  =  -  /  g  ^  dx  V  ^  H,^  H  J  . 

L  ^  i  6\i  i.  .'J  0 

t?t.  t  r,  f’ 

h  n 

Using  the  regularity  result  (4.2)  one  can  easily  show  that  if  .  is  the  solution  of  (4.1) 

and  u  r  -Fill  ,  then  (u,'i)  is  a  solution  of  (4.4),  and  if  (u,l)  is  a  solution  of  (4.4)  , 

then  i)i  is  the  solution  of  (4.1)  and  u  =  -iO  .  (4.4)  is  an  example  of  problem  ?  ir. 

section  2  with  (/^^  -  II  *11  -  ll-IlQ^h  '  “"h  "  "h  "o  '  " ‘V  “  "‘"z.h  '  “ 

/  u  V  dx,  and 
0 

b.  (u,sf)  «  y  /  u4  dx  -  /  u(J  )  ds  , 

n  -  ¥  m  t" 

TeT.  T  r. 

h  h 

(and  with  g  replaced  by  -g) .  Here  the  sulxrlass  of  data  for  which  (4.4)  is  uniquely 
solvable  is  D  »  0  *  H  ^(0)  . 

As  pointed  out  above,  can  be  identifies  with  (0)  ^  L^IT^) .  Under  this 

identification,  H^(n)  is  considered  a  linear  manifold  in  through  the  mapping 

n 

H^(n)  >  u  -  (u.uL  )  ^  L.fO)  •  L,(r.)  -  H°  . 

r.  ^  2  n  n 

ri 

Thus  an  element  u  ■  (u,u)  e  L.iFl)  •  L,(r.)  is  considered  to  be  in  H  (f!)  if  u  ^  H  (') 

2  2  n 

and  u|j,  »  u  .  To  be  completely  precise  bj^  should  be  defined  by 
h 

b.  (u,i#)  ■  y  /  uA^dx  -  /  u  (J  1^)  ds 

h  _  T  i  r.  »V 

Teth  T  Fj, 

for  u  »  (u,u)  e  “  L.,(n)  •  L-(r.  )  and  >f  t  }i^  .  Note  that 
n  2  2  n  n 

(4.5)  b.  ■  -/  7u  •  dx 

1  2  “ 

for  u  f  H  (n)  and  iP  e  H.  .  We  further  note  that  it  is  immediate  that  (2.1)  and  (2.2) 

n 

are  satisfied  with  constants  that  do  not  depend  on  h  . 

For  finite  dimensional  spaces  we  choose  ® 

is  defined  in  (3.3).  Problem  P.  thus  has  the,  form: 

n 
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Given  g  «  find  ‘  '-'j,  ^  '■’},  satisfying 


(4.6)  \ 


I  /  u  V  dx  +  7 

i  n  ^  Tff 


r  Viy^dx  -  f  v(J  -:7^)dx  =  0 
*  n  _  •  o  w 


)  r  u,  A  ;dx  -  I  u  (J  ■; — )ds  =  -  j  g  g  dx 
L  >  h  n  iv  i 

el,  T  .. 


Using  (4.5)  one  easily  sees  that  the  approximation  procedure  determined  by  (4.6)  is  tr.e 
same  as  that  considered  by  Glowins)ci  [14]  and  Mercier  [19]  and  further  developed  by 
Clarlet-Raviart  (9) .  Note  that  this  method  yields  direct  approximations  to  ,  and  to 
u  «  -Ai()  (the  streeun  function  and  vorticity  in  hydrodun2unical  problems)  . 

We  have  already  observed  that  (2.1)  and  (2.2)  are  satisfied,  in  order  to  apply 
Theorem  1  we  must  chec)c  the  stability  condition  (2.5)-(2.6). 

Theorem  2.  There  is  a  constant  >  0  ,  independent  of  h  ,  such  that 


1/  u  vdx| 


>  Y-  Hull 


V  u  «  z. 


i.e.,  (2.5)  is  satisfied. 

Proof.  Using  Lennia  1  we  have 

I  /  u  v  dx 

"'^Ml  V 

ve Z  0,h 


^  sup 


J  u  V  dx  I 

n 


0,h  Cllvllp 


=  c‘^  Hull, 


C"^  Hull, 


V  u  £  Z^  , 
n 


where  C  is  independent  of  h  .  Thus  (2.5)  holds  with  =  C 

Now  we  consider  (2.6).  Let  S,  =  {v  £  S.  :  f  v  dx  «  0). 

h  h  ^ 

Lemma  5.  There  is  a  constant  Cj^  >  0  ,  independent  of  h  ,  such  that 

1/  Vv  ’  V  dx  I 
o 

inf  sup  V  h  . 

VCS^  “'^"o,h"''"2.h 
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Proof.  We  first  note  that 


This  is  a  consequence  of  the  fact  that  an  operator  and  its  adjoint  have  equal  norms 


Given  vs  we  choose 


V  and  using  Lemma  1  we  obtain 


IS  convex 


is  the  Neumann  projection  of 


into  S  and  it  is  well-known  that 


0  on  r  we  see  from  the  definition  of  II  •!! 


from  (4.9)  that 


From  Lemma  4  with  k  =  1  and  r  «  2,  and  (4.10)  we  have 


sing  Lemma  2,  (4.9),  (4.10),  and  standard  approximabllity  results  we  find  that 


(•J.13) 


_<  C  h  (llv'  -  v'll^  +  lls:  -  ^11^) 


<  C  h  ^(hllv'll^  +  hlUllj) 


<  Cllvll^  . 


Now,  using  (4. 11)  -  (4 . 13)  we  )iave 
(4.14) 


<  C3  HvIIq  <  C3  llvllp^^ 


where  C3  is  independent  of  h  . 


Combining  (4.8)  and  (4.14)  we  get 

1/  Vv  .  dx| 


(4.15) 


inf  sup 


llvik  ^  llv'll. 


=  C3  >  0 


"h  "  h 

The  desired  result  now  follows  from  (4.7)  and  (4.15) . 

Theorem  3.  There  is  a  constant  )c„>0,  independent  of  h  ,  such  that 


(4.16) 


'b  (v,<)  I 

T  “MiTT- V''-\h 

0,h 

n 


V  e  W.  and  V  h  , 

n 


i.e.,  (2.6)  is  satisfied. 

Proof.  Let  V  f  and  set  e  =  j  ^  dx  .  Then  |e|  Clliifll^  and  if  =  f-e  e  S^. 


By  Lemma  5  there  is  a  v,  «  S.  such  that 

I  n 


bj^(Vi,.^)  =  bj^(V3,<)  =  -  /  VV3  •  Vif  dx  >_  Il.?ll2  ^  '''^''2,h  "  '■4"'^"o 


(4.17) 
and 

(4.18)  "Vo,h-'^''-^"2,h^  V'^"2.h 
We  also  know  that 


(4.19) 
and 

(4.20) 


-b^(v',y')  -  /  |7  >;rdx  >  Cgll^llg 


"^'"o.h^S  "^''2,h  • 
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Then,  using  (4 .17)- (4 .20) ,  we  have 


Combining  (4.21)  and  (4.22)  we  have  (4.16)  with 


We  are  now  ready  to  apply  Theorem  1  to  analyze  the  Ciarlet-Raviart  method.  We 


obtain 


r  ^  3,  and  suppose  1c  >  2.  Using  Lemmas  3  and  4  we 


In  addition,  (4.23)  yields  the  estimates 


We  now  derive  an  estimate  for 


by  means  of  the  well-known  duality  argu 


ment.  Given  d  e  H  (())  »  let  e  be  the  solution  of 


■A0  then  from  (4.2)  we  have 


Also,  from  the  discussion  following  equations  (4.4)  we  know  that  the  pair  (w,  ^) 


using  the  exact  equations (4 . 4) ,  and  the  -itz- 


Settinq  v  »  u  -  and  , 

Galerkin  equations  (4.6)  we  get 


(d,  0 


u,  ,w)  -  h.  (w.C’  -  )  -  b.  (u  -  u, 

h  h  h  h  h 

’1.  ,w  -  z)  -  b.  (w  ~  z,  .  -  .  )  -  b  (u  -  u, 

h  n  n  n 

w  (z,_)  'V 


Thus,  using  (2.1),  (2.2),  (4.26), 
l(d,^  -  <  Cdiu  - 

(4.27) 


and  Lemma  3  and  4  we  get 

+  II  ij)  -  'i.  II ,  v.)  llw-zl! 

h  2,h  , 


0,h 


+ 


II  u 


<  C  h  lldll  (llu  -  u.ll  ,  +  II 0  -  'v,ll-  ,  ). 

—-In  0,h  h  2,h 

Finally,  combining  (4.23)  and  (4.27),  we  have 

l(d,0-'i)! 

(4.28)  II  il;  -  =  sup  - -  1  C  h  ll  ■:■ll^ 

dfH~^(P.) 

where  s  =  min(r,)i+l). 

Estimates  (4.24)  and  (4.28)  improve  on  those  in  Ciarlet-Raviart  [9].  Scholz  [231 

obtained  (4.24)  under  the  assumption  that  ?  is  smooth.  (4.24)  and  (4.28)  were  also 

obtained  by  Fal)c-Osborn  [12].  Note  that  the  approach  of  this  paper  does  not  yield  error 

estimates  for  the  case  lc=l  for  the  method  studied  in  this  subsection  (and  also  for  the 

method  of  Subsection  b) ;  for  this  case  the  reader  is  referred  to  Scholz  [24].  Using 

L  -estimate  techniques  Scholz  [24]  has  shown  that  II  u-u,  II  „  =  0(h^  ^^^)  under  different 
*  ho 

assumptions  than  those  made  in  (4.24).  In  [25]  it  is  shown  that  in  any  subdomain 

Jl  cc  JJ,  IIu-ulII.  is  of  "nearly"  the  same  order  as  II  II  „  provided  .  is  suffi- 
U  n  h  U,i; 

ciently  smooth.  Finally  we  note  that  our  approach  allows  the  treatment  of  the  case  when 

2  -1 

g  €  For  example,  we  could  treat  the  case  where  g  is  the  Dirac  func¬ 

tion,  vrhich  corresponds  to  a  concentrated  load  in  plate  theory. 

Estimates  (4.25)  are  new  for  this  problem.  (4.25c)  provides  an  estimate  on  the  rate 
at  which  the  jumps  in  the  normal  derivatives  of  across  interelement  boundaries  is 

converging  to  zero  and  also  contains  the  estimate 


<  C  h 


s-3/2 


lUill^  . 


-21- 


b)  Hertnann-Miyoshi  method 

In  this  subsection  we  consider  another  mixed  method  for  the  approximate  solution  of 
(4.1).  In  this  method  the  auxiliary  variable  is  the  matrix  of  second  order  partial 
derivatives  of  .  . 

For  T  t  a,,(j  ^  *  (y  j  1  ^  j  <  j  ,  with  v,  .  •:  (T)  and  v  =  v_,  we 

h  ij  —  —  1]  12  21 

set 


and 


M  ( V)  * 


M  (v) 

VT  ' 


t2 

1. 

i.  j  =  l 


V.  .  V  .  > 
11  1  1 


v2 

L 

i.  j=l 


V  .  .  V  .  T  . 
11  1  1 


where  v  »  *'*l''^2'  unit  outward  normal  and  t  =  “  ^'^2'  ”^1^  unit 

tancjent  along  IT  .  We  note  that 

(4.29)  f  I  (V.  .  -  ^  +  — ii  )dx  =  /  (M  (v)  +  M  (v)  )ds 

i.j=l  T  ST  VT  -  3t 


for  all  c  f  H  (T) .  On 


('^(m  F  {V  =  (v..),  1  <  io  <2=  v^2  -  ...  . . 


=  V  ,  V. .  e  K  (T)  V  T  t  T  ,  and 
21  1  j  n 

II  (V)  is  continuous  across  interelement  boundaries} 


we  define 


llvll  .  =  y.  .  /  iv.  .  |“dx  +  h  /  |m(v)  I  ds, 

-  o,h  '•1.1  Q  11  r  ' 


1  2 

w^«re ,  on  an  interior  edge  T '  =  3T  n  3T  of  T  ,  we  set  M(v)  =  M  i (v)  =  M  o(v)  , 

n  -  v-*-  '  * 

and  on  a  boundary  edge  T*  of  T,  ,  we  set  M(v)  *  M  (v)  .  Then  we  define  (/. 

n  '  V  ^  h 

o 

to  be  the  completion  of  I'  with  respect  to  l!vl(  ,  .  It  is  clear  that 

h  '  0/h 

(1.30)  I'vll  .<{)■.  .  Ilv.  .11^ 

'  0,h  —  13  0,h 


for  all  V  <  H^(.')  =  {v  =  (v^j),  1  1  i.j  i.  ^  '^21'  '^ij 

norm  will  be  clear  from  the  context  whether  we  are  applying  it  to  scalar- 

2  1 

valued  or  matrix-valued  functions.  As  in  Subsection  a  vfe  let  ((',  =  H  n  H  .  Then  the 

h  h  0 

mixed  method  studied  in  this  subsection  is  Isased  on  the  following  formulation  of  (4.1): 
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iven  g  £  find  (u.i;,)  f  1/  x  ft'  satisfying 


Using  (4.29)  we  can  easily  establish  the  relations  between  (4.1)  and  (4.31).  If  is 


the  solution  of  (4.1)  and  u 


then 


(u,.)  is  a  solution  of  (4.31) 


is  the  solution  (4.1)  and  u 


then 


(4.31)  is  an  example  of  problem  P  with  I'  and  ft'  as  akiove 


Letting  S,  be  as  defined  in  (3.1),  we  consider  the  approximate  problem  P  with 


With  this  choice  for  the  forms  a  and  b.  and  spaces  V  and  W  ,  problem  P, 


describes  the  Herrmann-Miyoshi  method  [15,16,20  1.  Note  that  with  this  method  we  obtain 


(the  displacement  and  moments  in  elasticity 


direct  approximations  to  i  and 


problems) 


In  order  to  apply  Theorem  1  we  must  chec)c  (2.1),  (2.2),  (2.5),  and  (2.6).  (2.1) 


and  (2.2)  are  immediate.  In  light  of  (4.30),  the  proof  of  (2,5)  is  similar  to  the  proof 


of  (2.5)  for  the  method  in  Subsection  a.  Finally  we  consider  (2.6).  Let  (  W  be 


given.  By  Theorem  3  we  )cnow  there  is  a  v  £  S  such  that 


vie  immediately  nave  v  '  W, 


This  proves  (2.6). 

We  are  now  ready  to  apply  Theorem  1  to  analyze  the  Herrmann-Miyoshi  method.  This 
application  is  essentially  the  same  as  that  in  Subsection  a.  We  use  the  approximabilit; 
results  in  Lemma  3,  as  modified  for  matrix-valued  functions  with  the  aid  of  (4.30),  and 


in  Lemma  4.  We  will  just  state  the  results 


r  >  3.  Then 


where  s  =  min(r,)t+l).  From  (4.32)  we  obtain 


these  estimates.  We  further  note  that  (4.32)  contains  additional  information  corres 


pondiny  to  the  mesh  dependent  norms  (cf.  (4.25)) 


c) 


Hermann-johnson  method 


In  this  subsection  we  consider  a  further  method  for  the  approximate  solution  of  (4.1) 

in  which,  as  in  the  case  treated  in  Subsection  b  ,  the  auxiliary  variable  is  the  matrix 

of  second  order  partials  of  i(i  .  Also  as  in  Subsection  b  ,  the  method  is  based  on  the 

variational  formulation  (4.31)  (the  spaces  (/._  and  W.  and  the  forms  a,  and  b,  are 

h  n  n  n 

the  Scune  as  in  Subsection  b)  . 

We  now  consider  the  problem  with 


and 


Ww  *  s.  n  H^(Q) 
h  n  0 


C  P,,.!  V  T  e  ) 


This  choice  leads  to  the  Herrmann-Johnson  method  (15,16,17). 
differs  from  the  Hertmann^^iyoshi  method  only  in  the  choice  of 


space  V, 


h  ■ 


Note  that  this  method 
the  finite  dimensional 


This  example  has  certain  special  features  which  allow  an  analysis  that  is  rather 
different  than  that  employed  in  the  previous  two  examples.  These  special  features  involve 
the  existence  of  two  particular  projection  operators  denoted  by  and  .  We  turn 

to  this  now  . 

TT.  is  defined  as  in  [7,  Section  4).  For  v  =  (v.  .)  with  v,  .  e  H^(T)  and  v  = 
h  '  1]  1]  12 

v„,  we  define  'ir„v  »  (w..)  with  w.,  e  p,  ,  and  w, ,  =  w  by 
21  T-  i]  ij  k-1  12  21 


(4.35) 


/  M  (V  -  IT  v)f  ds  «  0  V  f  €  P,  ,  and  for  each  side  T'  of  T  , 
i,  V  '  T'  K-1 


/  dx  -  0  V  f  €  P^_ 


2  • 


By  Lemna  3  in  Ifi] ,  D^V  is  uniquely  determined  by  (4.35).  Now  for  v  «  we  define 
■h  ‘  "h 

'■>.-”1  ■  't'-l  ’ 
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since  we  can  write 


it  is  clear  that 


concerning  the  approximation  of  y  by  ■"«,  y  have 


r  >  3 .  Then 


it  is  shown  that 


Thus  it  remains  to  show  that 


and  assume  T  is  the  image  of  T  under  the  mapping  F(x)  =  Bx  +  b 


Given  a  matrix  valued  function  w(x)  on  T  we  set  w(x)  =  C  w  (F(x))  C 


C  =  B  .  (Note  that  the  correspondence  between  (matrix  valued)  functions  on  T  and  on  T 
is  different  than  the  one  introduced  in  Section  3.)  Recall  that  v  =  C^v|b*'v|  ((3.5)). 


Then  we  have 


max 


C(f)h1  lfl.4  21-1  ,2 

2(:+2)  ’’  ’-U.T 


This  completes  the  proof. 


The  second  projection  operator  is  the  interpolation  operator  1^^  introduced  in 

®  2  1 

Section  3.  As  in  the  proof  of  Lemma  5  in  (7]  ,  for  y  e  and  e  H  (f2)  n  we 

can  write  _ 

3  V .  . 

(4.38)  b.  (v,.p)=  -Y  J-  .  /  t — <3x  +  I  /a(T',v)i(>  is  I  B(a,v)^(a) 

"  '  T  T  3*i'**j  T'er  T*  '  aeJ^ 

n  n 

where  is  the  set  of  all  sides  of  the  triangulation  ,  J.  is  the  set  of  all  vertices 

h  h  n 

of  T,  ,  and  A(T',v)  is  a  polynomial  of  degree  less  than  or  equal  to  k-2  in  the 
h 

3V.  , 

variable  s  .  Since  for  v  e  Vj^  we  have  1  t  ^k-3  A(T',v)  e 

follows  from  (4.38)  that  Z,v>  =  I.  'f ,  as  defined  in  Section  3,  satisfies 

h  h 

(4.39)  bj^(v,  2:j^s<'-sJ)«0  VveVj^. 

Now  we  are  ready  to  derive  the  error  estimates .  First  we  estimate  Hu-  Uj^ll  ^  . 
Subtracting  (2.4a)  from  (2.3a)  we  obtain 

(4.40)  “h^-  ”  *  *’h*-  '  ”'*’h^  ”  °  ^  *  '^h  ■ 

Suppose  V  f  Zj^  s  {w  e  Vj^:bj^(w,^)  «  0  V  ip  e  Mj^).  Then,  from  (4.39)  we  see  tliat 

b..(v,«)  »  »  0  for  all  'fi  e  H^(fl)  o  H^CO).  Hence  from  (4.40)  we  have 

h  '  h  '  h  u 

(4.41)  a.  (u  -  u.  ,v)  »  0  V  v  e  Z.  . 

h  '  -h  '  '  h 

Subtracting  (2.4b)  from  (3.3b)  and  using  (4.36)  we  see  that 

*’h*’'h'  ”  ”  ^h*-  *  °  ^  ^  ’^h  ' 

i.e.,  "t.u  -  u.  f  Z.  .  Thus,  recalling  (4.41), 
n  -n  n 


yi 


11“  -  “wIL  “  ~  “v.»“  ■ 

^  -h  0  h  ,  .h  . 


a.  (u  -  u,  ,(u  -  u.  )  -  (tt.u  -  u  )) 
n  -  -h  "  “h  n'  n 


<  II U  -  U.  II  II U  -  11.  ull 
—  '  -ho-  h-  0 


and  hence 


II U  -  U.ll„  <  II  IT.  U  -  ull 
'  -h  0  —  h'  '0 


Suppose  now  that  £  H  (fl)  ,  r  3  .  Then  from  (4.42)  and  Lemma  6  we  have 


llu-^lloiCh=ll4,ll^^2 


where  s  »  min  ()c,r-2) . 


Now  we  estimate 


.  As  in  Subsection  a  we  can  write 


(4.44)  n  ”  -a.  (u-u.  ,w-z)  -  b.  (w-z,i^-\|i.  )  -  b  (u-u  ,6-u)  V  {z,\i)  c  V  « 

no  n-^n*--  n--  n  n-'n  -  n 


where  6  is  the  solution  of 


and  w 


ij  3x^3x^ 


(4.45) 


A  6  =  d  e  Lj  on  SI 

6  »  -1^  «  0  on  r 
3v 


We  note  that  (w,6)  satisfies 


a.  (w,v)  +  b,  (v,e)  »  0  V  V  £  (/. 
h  -  -  h  -  '  h 


b.  (w,^)  »-/  dv>dx  V  sfflll. 


(cf.  (4.31)).  In  (4.44)  let  2  =  it.  w  and  y  =  £  6  .  This  gives 

n  ”  n 

(4.46)  -  -aj^(u  -  V  ^  ■  V’  ■  ’=h''?'V''‘'"V  ■  ‘’h‘-"'Jh'®‘^h®'  • 


V7e  now  estimate  each  tern  in  (4.46). 


Using  (4.36),  (4.39),  (4.45),  and  Lemma  3  we  have 


(4.47) 


(b^(w,t)i-rj^i))) 


«  I  (d,i|(-£j^i)))^| 

<  lldllgll,^-EJ^tllg 

<  C  h®  lUII  II  dll 


where  b  «  nln(r-l,k4>l)  • 


^Hi 


In  our  estiAate  for  the  third  term  on  the  right  side  of  (4.46)  we  treat  the  cases 


First  assume  k  >  2.  Then,  using  (4.26)  ,  (4.39) ,  and 


LeMMs  4  and  6,  we  find  that 


where  s  «  min(r-l,k+l)  .  Mow  suppose  k  =  1.  Then.using  (4.26),  (4.31),  (4.39),  and 


Leiana  3  we  have 


Finally,  using  (4.26),  (4.42),  and  Lensna  6  we  obtain 


where  s 


Conbining  (4. 46) -(4. SO)  we  have 


One  can  also  prove  that 


Estimate  (4.53)  improves  on  estimates  in  (7],  Estimates  (4.43),  (4. 51)- (4. 53) 
are  given  in  (7),  and  (4.43)  and  (4.51)-(4.54)  are  proved  in  [12]. 

Remarks :  1)  As  in  Subsection  b  we  could  have  shown  that  the  method  studied  here  is 

stable  with  respect  to  the  norm  II  II  „  .  +  II  II  -  .  ,  and  then  obtained  error  estimates  in 

0  /  /  n 

this  norm.  This  approach  would  have  allowed  the  treatment  of  the  case  vhen 
2  -1 

g  €  (Hj^  )'  -  H  (n)  (cf.  the  next  to  the  last  paragraph  in  Subsection  a).  However, 

due  to  the  special  nature  of  this  example,  more  refined  estimates  can  be  obtained  by 
the  analysis  sketched  above  in  the  case  when  sufficient  regularity  of  the  solution  is 
assumed.  Thus  the  mesh  dependent  norms  play  a  less  central  role  in  the  ^malysis  of 
this  method  th^ul  in  previous  methods.  They  are,  however,  convenient;  their  use  leads 
to  a  very  natural  setting  for  the  study  of  this  example. 

2)  The  analysis  in  this  subsection  was  based  on  the  projections  7t^  and  2md 

the  fact  that 

2  c  2  S  {w  :  w  £  (/  ,  b.  (w,^)  ”0  V  ifi  €  H^(n)  n  H^(n)} 

h  -  '  h  h  -  0 

which  follows  from  the  existence  of  Z.  .  For  a  general  discussion  of  the  projections 

n 

ir.  £md  Z.  and  the  condition  z.  c  z  see  Falk  and  Osborn  [12]  and  Fortin  [13]. 
n  h  h 

3)  In  this  subsection  the  mesh  family  is  not  required  to  be  quasi-uniform. 
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